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Abstract 

A sailplane being developed at NASA Dryden 
Flight Research Center will support a high-altitude flight 
experiment. The experiment will measure the 
performance parameters of an airfoil at high altitudes 
(70,000 to 100,000 ft), low Reynolds numbers (200,000 
to 700,000), and high subsonic Mach numbers (0.5 and 
0.65). The airfoil section lift and drag are determined 
from pitot and static pressure measurements. The 
locations of the separation bubble, Tollmien-Schlichting 
boundary layer instability frequencies, and vortex 
shedding are measured from a hot-film strip. The details 
of the planned flight experiment are presented. Several 
predictions of the airfoil performance are also presented. 
Mark Drela from the Massachusetts Institute of 
Technology designed the APEX- 16 airfoil, using the 
MSES code. Two-dimensional Navier-Stokes analyses 
were performed by Mahidhar Tatineni and Xiaolin 
Zhong from the University of California, Los Angeles, 
and by the authors at NASA Dryden. 

Nomenclature 
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c chord, ft 

C d section drag coefficient 

Cj section lift coefficient 

C m section moment coefficient 

C p pressure coefficient 

du derivative of velocity, ft/sec 

dx derivative of length, ft 

e" amplitude ratio 

EMI electrical magnetic interference 

g force of gravity 

gm gram 

KEAS knots equivalent airspeed at sea level 

M Mach number 

n rit critical amplification parameter 

P pressure, lb/ft 1 2 

PCM pulse code modulation 

P max dimensionless pressure gradient 

P s static pressure, lb/ft 2 

Pj total pressure, lb/ft 2 

P.T. pressure transducer 

q dynamic pressure, lbm/ft/sec 

Re Reynolds number 

RF1 radio frequency interference 

RTD resistive temperature device 
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sep separation 

U uncertainty of the variable x 

x TR transition location 

a angle of attack, deg 

y ratio of specific heats 

0 momentum thickness 

u kinematic viscosity, ft/sec 2 

Introduction ^- 

The need for cost-effective high-altitude vehicles to 
conduct atmospheric research has created interest in 
high-altitude (low-Reynolds-number) airfoils. In 
support of this need, NASA Dryden Flight Research 
Center is developing a sailplane called APEX that will 
measure the parameters affecting the performance of the 
airfoil in actual high-altitude flight. The APEX sailplane 
will be released from a high-altitude balloon from 
approximately 108,000 ft altitude and then remotely 
piloted. Figure 1 shows a schematic of the flight 
mission. 



'iso of trade names or names of manufacturers in this document 
does not constitute an official endorsement of such products or 
manufacturers, either expressed or implied, by the National 
Aeronautics and Space Administration. 


The first 30 sec after release from the balloon are the 
most critical for the APEX flight control system. 
Transition to horizontal flight occurs during this period 
with the assistance of four small rockets, which have a 
combined thrust of 784 lb. After the transition to 
horizontal flight, the airfoil parameters affecting 
performance are measured as the sailplane descends 
from 100,000 to 70,000 ft. The sailplane is then brought 
to a horizontal landing at the Rogers dry lakebed at 
Edwards Air Force Base, California. 

Low-Reynolds-number airfoils typically exhibit 
laminar separation bubbles as shown schematically in 
figure 2. These separation bubbles are known to 
significantly affect the performance of an airfoil. The 
bubble is formed when the laminar flow separates as a 
result of encountering the adverse pressure region of the 
airfoil. The separated free shear layer is unstable, which 
amplifies the Tollmien-Schlichting instability waves. 
The free shear flow generally transitions rapidly from 
laminar flow to turbulent flow and then reattaches to the 
airfoil surface. The lambda shocks, which occur in the 
transonic flight regime, are expected to increase the 
amplification of the Tollmien-Schlichting instability 
waves. 

The objectives of the APEX experiment are 

• To increase the understanding of airfoil 
performance in the high-altitude, low-Reynolds- 
number, and high-subsonic-Mach-number flight 
regime. 

• To obtain flight test data of airfoil performance 
parameters that can be used for validation of airfoil 
design codes. 



gradient 


Figure 2. Laminar separation bubble. 
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This paper presents a description of the APEX 
experiment. The design details used to determine the 
perfonnance parameters of the airfoil will be discussed. 
A preliminary error analysis will also be presented. 
Various numerical models used to predict the airfoil 
performance parameters will also be discussed. 

Previous Research 

Several studies investigating the perfonnance and 
characteristics of low-Reynolds-number airfoils have 
been performed. Mueller 1 presents an excellent 
summary of the research before 1985. One interesting 
aspect that Mueller discusses is the hysteresis that often 
occurs in the drag polars. Mueller’s wind-tunnel studies 
show that the airfoil perfonnance, including the 
hysteresis, could be significantly affected by free-stream 
turbulence and surface roughness. LeBlanc et al. 
performed wind-tunnel measurements on a Liebeck 
airfoil at low Reynolds numbers. The boundary layer 
velocity fluctuations in the separation bubble were 
measured with hot-wire anemometry. The measured 
velocity spectra of peak instability frequencies agreed 
with tlie predictions from linear stability theory 
analysis. 

Howard and Kindelspire 3 performed a wind-tunnel 
study of the free-stream- turbulence effects on an airfoil. 
Their investigation shows that transition develops more 
rapidly by increased free-stream turbulence and also 
that boundary layer instability growth is greatest when 
the turbulence length scale is on the same order as the 
boundary layer thickness. Dovgal et al. 4 discuss several 
aspects of the instability associated with laminar 
separation bubbles including receptivity, linear 
instability, and nonlinear interactions. 

Recent investigations of Pauley et al., 5 Ripley and 
Pauley, 6 and Muti Lin and Pauley 7 show that the 
separation bubble may become unstable under certain 
conditions, and experiences periodic vortex shedding. 
Their transient incompressible Navier-Stokes analytical 
studies characterized the unsteady vortex shedding 
structure. Tatineni and Zhong 8, 9 performed a two- 
dimensional, time-accurate Navier-Stokes analysis on 
the APEX- 16 airfoil flow field. Their analysis indicates 
that the airfoil separation bubble is unstable and 
periodically sheds at the flight conditions of the APEX 
sailplane. A linear stability analysis also showed that the 
most dominant instability frequency matches the 
frequency of the periodic vortex shedding. 

Drela 10 investigated high-altitude, low-Reynolds- 
number airfoils in the transonic flight regime with the 


ISES code. An interesting aspect of this investigation is 
that airfoil performance in the high-altitude flight 
regime may depend largely upon the effectiveness of 
lambda shocks to increase the amplification of 
instability waves and increase the transition rate in the 
separation bubble. A conclusion is that experimental test 
data exploring laminar shock-boundary layer 
interactions and its effect on Tollmien-Schlichting 
instability waves would significantly reduce the 
uncertainties in the analysis. 

Turbulence is expected to be a major factor in the 
performance of the APEX- 16 airfoil. The length of the 
separation bubble depends on the growth of the 
instability waves within the free shear layer and 
transition to fully developed turbulence. The process by 
which free-stream turbulence enters the boundary layer 
and becomes amplified is known as receptivity. 
Qualitatively, the concept is simple: Free-stream 
turbulence of various amplitudes and wavelengths enters 
the boundary layer and either decays or grows. The most 
unstable wavelengths grow at the greatest rates and 
develop into Tollmien-Schlichting instability waves and 
eventually cause transition to turbulence. Generally, the 
larger the free-stream turbulence amplitudes, the 
quicker the boundary layer transitions to turbulence. 
Quantitatively, receptivity is a complicated subject that 
is not well understood. Reshotko" summarizes the 
current understanding of receptivity. Although the 
receptivity process to transition is difficult to 
quantitatively predict, it is well known to be a strong 
function of free-stream turbulence as shown by 
Dry den et al. 12, 13 in several experiments measuring the 
critical Reynolds number of a sphere as a function of 
free-stream turbulent intensity. Fisher and Dougherty 14 
performed a series of transition measurements on a cone 
in wind tunnels and in flight. Then results show that the 
transition location is a function of the free-stream 
turbulence. The atmospheric turbulence in flight is very 
low and of large wavelength in relation to the thickness 
of the boundary layer. 

Presently, no existing wind tunnel can provide the 
high-altitude (70,000 to 100,000 ft), low-Reynolds- 
number (200,000 to 700,000), high-subsonic-Mach- 
number (0.5 to 0.65), and low-free-stream-turbulent- 
intensity (0.02 percent or less) environment necessary to 
accurately measure the APEX- 16 airfoil perfonnance. 
Natural atmospheric turbulence is the rationale for 
constructing the APEX research sailplane and 
measuring the in-flight airfoil performance parameters 
rather than performing a wind-tunnel study. 
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APEX Sailplane Description 


Table 1. Coordinates (x, y), inches. 


Murray et al. 15 originally proposed the APEX 
sailplane as a modified Schweizer SGS 1-36 sailplane. 
As the APEX design proceeded, the sailplane evolved 
into the current configuration. Figure 3 shows the APEX 
sailplane geometry. The sailplane is 22.7 ft long with a 
wingspan of 41.2 ft and has a whig aspect ratio of 13.6. 
The experimental test section, where the performance 
parameters ate measured, is at the midspan point of the 
right wing as shown in the figure. The sailplane is 
designed for a target gross weight of 600 lb with a 5 -g 
maneuver load factor. The airframe has been designed 
and is currently being fabricated by Advanced Soaring 
Concepts (Camarillo, California) from graphite/epoxy 
and boron/epoxy composites. 

Drela designed the APEX- 16 airfoil using the coupled 
viscous/inviscid MSES code. 16 ' 17 The coordinates and 
airfoil shape are shown in table 1 and figure 4, 
respectively. The airfoil dimensional tolerances for the 
wing construction are specified at ±0.005 in. to reduce 
the effects of surface roughness and waviness. The wing 
is a rectangular planform of the APEX- 16 airfoil, as was 
shown in figure 3. The whig incorporates a 2-deg linear 
washin to reduce three-dimensional (spanwise) effects 
and to provide a more uniform section lift coefficient 
( C i ) distribution over the experimental area of the whig. 


Area 


,101.21 in., 


a 




77.70 74.40 
in. in. 


Wing 

124.61 ft' 

Aileron 

9.52 ft' 

Horizontal 

stabilizer 

20.35 ft ; 

Elevator 

6.80 ft' 

Vertical fin 

11.87 ft' 

Rudder 

3.62 ft' 





airfoil 1 


1.5° 

61.18 in. \ 

78.03 in. 
-< = 

■ 

i 
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Figure 3. APEX three-dimensional view. 
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Figure 4. APEX-16 sectional (37. 22-in. chord). 


The forward undercut camber on the lower surface does 
not directly affect the separation bubble and was 
incorporated into the design to increase the maximum 
lift coefficient and decrease the pitcliing moment. The 
airfoil was designed to provide good performance 
throughout the entire APEX flight envelope. Figure 5 
shows the predicted APEX flight envelope. The 
challenge of the design was to correctly predict the 
characteristics of the separation bubble without 
experimental data for code verification hi the high- 
subsonic-Mach-number and low-Reynolds-number 
flight regime. 

The experiment is limited by several design 
constraints. Weight is a major design consideration for 
the experiment. The gross vehicle flight weight for the 
sailplane is specified at 600 lb to provide an adequate 
stall margin for attaining a ceiling altitude of 100,000 ft. 
The experiment is limited to 10 percent of the gross 
vehicle weight or 60 lb. Packaging the experiment is 
also another major design constraint. Much of the 
instrumentation electronics, including pressure 
transducers, accelerometers, and the hot-film 
anemometry, signal conditioning, and analog-to-digital 
(A/D) conversion cards, are in the whig next to the 
experiment. The electronic instrumentation is located 
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Mach number 

990004 

Figure 5. APEX flight envelope. 

close to the experiment to minimize noise from radio 
frequency interference (RFT) and electrical magnetic 
interference (EMI). The whig chord is 37.22 in. with a 
maximum ulterior thickness of 5 in., which makes 
instrumentation packaging difficult. Another design 
constraint is the high altitude. The available off-the- 
shelf instrumentation that can provide adequate range 
and accuracy at high altitude is limited. Some 
instrumentation had to be specifically designed for the 
experiment. The air density at 100,000 ft altitude is 
approximately 1 percent of its value at sea level, which 
substantially lowers the convective cooling rates of the 
electronics to the air. Some electronics require a 
specialized cooling design to avoid overheating. 

APEX Experimental Description 

The flight experiment to measure the performance 
parameters of the APEX- 16 airfoil consists of three 
primary measurements: 

• First, to measure section lift a series of static 
pressure taps circle the airfoil at one spanwise 
location. 

• Second, to measure section drag a trailing rake sits 
behind the airfoil with a support sting. 

• Third, to measure the separation bubble location, 
Tollmien-Schlichting frequencies, and vortex 
shedding a hot-film strip sits on the top surface of 
the airfoil. 

Figure 6 presents a schematic layout of the airfoil 
instrumentation. In addition to the primaiy 
measurements, the instrumentation also includes a Kiel 
probe to measure free-stream total pressure, a trailing 


static probe to measure free-stream static pressure, a 
boundary layer rake to determine the velocity profile, a 
total temperature measurement, five integrating 
boundary layer rakes to determine the section drag 
developing over the upper surface, two integrating 
trailing rakes to determine total section drag, two vanes 
to measure angles of attack and sideslip, and three 
accelerometers to measure whig surface vibration. 


( 



Figure 6. APEX instrumentation layout on right-whig 
test section. 


Pressure Measurement System 

Figure 7 presents a schematic of the pneumatic 
pressure measurement system. Fifty static pressure ports 
along the chord (30 on the upper surface and 20 on the 
lower surface) measure the pressure distribution over 
the airfoil. The ports have a 0. 05-in. diameter and are 
staggered at a 15 -deg angle relative to the chord to 
prevent contamination from upstream ports. A trailing 
rake comprises 26 total pressure probes and 3 static 
probes to determine section drag ( C fI ). The rake is 
mounted 0.3 chord length aft of the airfoil where the 
static pressure is expected to be fully recovered. 

The airfoil section drag is calculated from the rake 

i o 

pressures based upon the Jones method corrected for 
compressibility effects. A Kiel probe, located mid- 
chord, 8 in. from the lower surface of the airfoil, 
measures a reference total pressure (P T ). A trailing 
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Figure 7. Pneumatic pressure measurement system for 
the right-wing test section. 


chosen because it is flight proven, lightweight (3.7 oz), 
and small (1 by 4 by 1/2 in.). Each transducer has 
32 input pressure port channels. The input pressure is 
differentially measured against a reference pressure. 
The transducer may have only one reference pressure, 
either the static pressure from the trailing static probe or 
the total pressure from the Kiel probe as was shown in 
figure 7. 

The absolute pressure is measured by two Baratrons 
transducers manufactured by MKS Instruments 
Incorporated (Andover, Massachusetts). The transducers 
measure the static pressure in two ranges, 0 to 28 lb/ft 2 
and 0 to 280 lb/ft 2 , and provide an absolute reference 
pressure measurement for the entire system. Both units 
are saturated below a 50,000-ft altitude, which is 
acceptable because the experiment is designed for high 
altitude between 70,000 and 100,000 ft. The data from 
all transducers are sampled and sent by telemetry to the 
ground at a rate of 25 Hz. The pressure system mounted 
in and on the right wing including rakes, probes, and 
transducers weighs less than 30 lb. 


static probe, placed three chord lengths aft of the airfoil, 
measures a reference static pressure ( P ). A 
conventional boundary layer rake placed at 70-percent 
chord determines the boundary layer velocity profile 
development on the upper surface. Five integrating 
boundary layer rakes are located at 60-, 70-, 80-, 90-, 
and 100-percent chord, and two integrating trailing 
rakes are located at 0.3 chord lengths aft of the airfoil. 
The integrating rakes are multi-pitot probe rakes such as 
developed by Silverstein and Katzoff 19 hi which the 
pitot probes are plumbed into a common reservoir for a 
single average total pressure measurement. The average 
total pressure measurement has been shown by Drela # to 
be a direct determination of the sum of momentum 
thickness and displacement thickness from which drag 
can be calculated. 

The pressure measurement system design uses four 
electronically scanned differential pressure transducers 
manufactured by Pressure Systems Incorporated 
(Hampton, Virginia) to measure differential pressure. 
The transducers are silicon peizoresistive pressure 
sensors with a range of ±52 lb/ft 2 . This transducer was 


Uncertainty Analysis 


A preliminary measurement uncertainty analysis was 
performed on the pressure system and is summarized 
with a discussion of the bias error for the calculation of 
lift. This analysis is based on the general uncertainty 
analysis of Coleman and Steele. 20 The pressure 
coefficient ( C ) can be defined in terms of measured 
quantities as follows: 


C 


P 


P~P S = P~P S 

4 0.7P S M 2 


P-P 


S 


0.7 P, 


y- 1 


Y - 1 -, 



( 1 ) 


where P is the port pressure, q is the dynamic pressure, 
M is the Mach number, and y is the ratio of specific 
heats. 


Treating P - P s , P , and Pj/P s as measured 
# Drela, Mark, “Integrating Rake Design,” private communication, quantities, the general uncertainty equation is 
1995 . 
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where U is the uncertainty for the subscript variable. 

For a C = 1 and a flight condition of Mach 0.65 
and 100,000 ft altitude, the expected pressure quantities 
and their associated bias errors from the manufacturer’s 
specifications are 

P-P s = 7.49 ±0.14 lb/ft 2 

P s = 22.7 ± 0.072 lb/ft 2 (3) 

P T /P S = 1.33 ±0.0012 lb/ft 2 

Substituting these quantities into the uncertainty 
equation (2) yields the pressure coefficient bias error of 


C p = 1 ±0.02 (or ±2%) (4) 

For a C p = -0.5, the result is 

C p = — 0.5 ± 0.02 (or ± 4%) (5) 

The bias error for the section lift coefficient, 

C 

c r;K-y* < 6 > 

0 

is expected to be 

C, = 0.9 ±0.028 (or ±3%) (7) 


where c is the airfoil chord and the subscripts / and u 
refer to the lower and upper surface, respectively. 

A similar bias error analysis for the calculation of the 
section drag coefficient for the flight condition of 
Mach 0.65 and 100,000-ft altitude yields 

C d = 0.02 ± 0.004 (or ± 20%) (8) 


The percentage of bias error decreases as the sailplane 
descends to lower altitudes because the measured 
pressures become larger. For a C p = 1 and a flight 
condition of Mach 0.65 and 70,000-ft altitude, the 
expected bias errors are 

C p = 1± 0.005 (or ±0.5%) (4a) 

C, = 0.9 ±0.006 (or ±0.7%) (7a) 

C d = 0.01 ±0.0005 (or ±0.5%) (8a) 

Marchman 21 states that “[djrag has always been the 
most difficult aerodynamic force to measure and the low 
drag forces occurring in low Reynolds number flows 
make the problem even more difficult.” The APEX 
experiment is not immune to these difficulties. The 
average wake deficit pressure is estimated to be on the 
order of 1 lb/ft 2 (M = 0.65, altitude = 100,000 ft) and 
5 lb/ft 2 (M = 0.65, altitude = 70,000 ft). The range of 
the differential pressure transducer used to measure the 
wake deficit pressure is ±52 lb/ft 2 and has an accuracy 
of ±0. 1 percent full scale. A smaller range transducer to 
lower die bias percentage error was not available that 
would satisfy the APEX design constraints of size and 
weight. The transducer is capable of remote zero-point 
calibration by applying a single reference pressure to 
both sides of the differential. This zero-point calibration 
will be performed before release from the balloon at 
100,000 ft altitude and just after the experimental 
measurements are completed at 70,000 ft altitude. 

A comparison was made by Marchman of file drag 
coefficient measurements on a Wortmann FX63-137 
airfoil performed by three different research facilities. 
The results show differences of more than 50 percent in 
the measurement of drag coefficient for similar test 
conditions. Although the section drag coefficient bias 
error range of ±20 to ±5 percent is large, it is not 
unreasonable given the difficulties of the measurement. 
The 5-percent error at 70,000 ft (Re = 700,000) is more 
representative of measured in-flight drag errors at high 
Reynolds numbers. For example Arnaiz 22 measured the 
in-flight drag on the XB-70 airplane accurate to 
±6.5 percent at Reynolds numbers ranging from 
1,000,000 to 3,000,000 per foot. 

Other errors are associated with the accuracy of the 
pressure system measurements. The displacement 
effects associated with the pitot probe disturbing the 
local flow field can affect the accuracy of the probe total 
pressure measurement. Montoya et al. discuss these 
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displacement effects in detail. Allen 24 reports a survey 
of pitot probe displacement corrections for free shear' 
flows of approximately 18 percent of the probe 
diameter. The APEX trailing rake pitot probes have a 
0.07-in. diameter and the airfoil wake height is expected 
to be 3 to 5 in. depending on flight conditions. The 
displacement effects of the pitot probe total pressure 
measurements are less than 0.2 percent. 

Lag times in the pneumatic tubing between the 
pressure port and the transducer can also be a large 
source of error under dynamic conditions. The low static 
pressure at 100,000 ft altitude, 23 lb/ft 2 , can cause large 
pneumatic lag times. Laboratory tests and calculations 
performed at NASA Dryden show that the lag times are 
less than 1.5 sec. Lag times are reduced as the sailplane 
descends in altitude because lag times are inversely 
proportional to pressure. Lag times should not pose a 
significant source of error as flight simulations show 
that a flight condition of constant Mach number and 
constant lift coefficient for the sailplane can be 
maintained for 5 to 10 sec. 

It is difficult to estimate how much pressure 
measurement uncertainty is a result of the three- 
dimensional effects of turbulence. Large-scale turbulent 
vortical structures and vortex shedding, if present, may 
significantly affect the accuracy of the static and pitot 
pressure measurements. The pneumatic tubing lags 
dampen the dynamic response of the pressure 
measurements. Pressure variations caused by the large 
vortical structures moving across the airfoil flow field 
are not detected and their effects on a time averaged or 
lagged static and pitot pressure measurements are 
uncertain. 

Large-scale vortical structures also present 
uncertainties caused by changes in the local flow 
direction. As shown by Mueller, 1 the accuracy of 
trailing rakes to determine drag is severely 
compromised in flows with large-scale vortical 
structures resulting from the changing flow direction. 
The changing flow direction may also affect a static 
pressure measurement in similar fashion as a vortex 
passes over a static port on the surface of the airfoil. The 
presence of large-scale turbulent vortical structures and 
vortex shedding is detected from the hot-film 
measurements. 

The experiments of Batill and Mueller 25 and 
Guglielmo and Selig 26 suggest that spanwise variations 
may exist hi the separation bubble dynamics caused by 
the three-dimensional effects of turbulence. Bastedo and 
Mueller measured a significant spanwise variation in a 


laminar separation bubble due to wing tip vortices. The 
APEX wing incorporates a 2-deg linear washin to 
reduce spanwise effects of the wing tip vortex and to 
provide a more uniform Cj distribution over the 
experimental area of the wing. 

Random errors are often difficult to separate from the 
natural in-flight pressure fluctuations. Random errors 
are believed to be largely caused by RL1 and EMI. 
Natural pressure fluctuations are caused by atmospheric 
turbulence, aircraft vibration, and changes to flight 
conditions such as angle of attack and angle of sideslip. 
Before releasing the sailplane from the high-altitude 
balloon, a series of pressure data samples is taken and 
sent by telemetry to the ground. An assumption is that, 
while the sailplane is suspended under the balloon, the 
air relative to the sailplane is still and there are no 
natural pressure fluctuations. All fluctuations in the 
prelaunch data samples, therefore, are assumed to be 
random errors. Any additional random fluctuations that 
occur in flight are assumed to be natural pressure 
fluctuations occurring over the surface of the airfoil. In 
addition, the pneumatic tubing lags dampen all high- 
frequency natural pressure fluctuations. 

The difficulty in obtaining accurate pressure 
measurements is realized when considering all the 
possible errors (bias, nonlinear pitot effects, pneumatic 
lags, large-scale vortical structures, spanwise variations, 
and random). The presence of large-scale vortical 
structures severely compromises the measurement 
accuracy. However, identifying flight regimes that 
contain large-scale vortical structures is an important 
pail of the experiment. As will be discussed later, these 
regimes are highly undesirable because of their 
associated large drag. The APEX experiment is 
expected to significantly increase the understanding of 
low-Reynolds-number airfoils at high altitude and 
provide data for validation of airfoil design codes. The 
uncertainties in the experimental data will be properly 
accounted for. 

Hot-Lilm Measurement System 

A multi-element hot-film strip is mounted over the 
APEX- 16 airfoil. The hot-film strip measures the state 
of the boundary layer (i.e., laminar boundary layer, 
laminar separation, bubble region, turbulent 
reattachment, turbulent boundary layer, turbulent 
separation, and vortex shedding) and the frequency of 
the Tollmien-Schlichting instability waves in the 
separation bubble. The strip consists of 50 hot films on 
the top surface in 2-percent chord increments starting at 
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zero-percent chord. The hot films ate spaced 
approximately 0.75 in. apart. In addition, four hot films 
are placed on the bottom surface at 10-, 30-, 60-, and 
90-percent chord. This hot-film-strip configuration is 
used for the initial flights. After the separation bubble is 
located for the APEX flight regime, the strip is replaced 
with a denser strip concentrated in the area of the 
separation bubble (approximately 50 evenly spaced hot 
films on a 15-percent chord length). 

The desire to obtain valid hot-film anemometry data 
at frequencies up to 10 kHz combined with the weight 
and packaging limitations significantly influenced the 
design. The APEX telemetry system cannot handle these 
high data transfer rates for the large number of channels. 
Therefore, the data are stored on board in random access 
memory and later, after the high-altitude tests are 
completed, sent by telemetry to the ground at lower data 
transfer rates. The hot-film data are split into two 
components: a DC component and an AC component. 
The DC component is sampled at 200 Hz and sent by 
telemetry to the ground in real time. The AC component 
is sampled in 1-sec data intervals at 20 to 25 kHz and 
stored in memory. The system is capable of storing up to 
ten 1-sec data intervals during a flight. The system is 
commanded from the ground uplink to begin storing a 
1-sec data interval of AC data. In summary, all 54 hot 
films are sampled for their DC component at 200 Hz as 
the sailplane descends from 100,000 to 70,000 ft. The 
54 hot films are sampled for their AC component at 20 
to 25 kHz in 1-sec data intervals for up to 10 intervals as 
the sailplane descends from 100,000 to 70,000 ft. 

Figure 8 shows a schematic of the hot-film 
anemometry and data acquisition system. The 
anemometry is a temperature-compensated system 
developed by Chiles. The anemometry incorporates a 
temperature sensor located in proximity to the hot film 
as part of the Wheatstone bridge. The temperature 
sensor corrects the sensitivity of the anemometry circuit 
to changes hi the average adiabatic wall temperature 
caused by changes in flight condition. The system also 
incorporates automatic gam ranging that automatically 
adjusts the AC signal gam so that the peak-to-peak 
signal is approximately 50 percent of the total signal 
range. Each hot film has an anemometry circuit, a signal 
conditioning circuit that includes the automatic gain 
ranging, and a 10-kHz low-pass filter to prevent anti- 
aliasing. 

The hot-film signals are multiplexed in groups of four, 
to a data logger and A/D converter (fig. 8). The data are 
then transferred to the pulse code modulation (PCM) 
encoder at lower data rates to be sent by telemetry to the 


Hot films (qty 54) 
(50 on upper surface) 



Figure 8. Hot-film anemometiy system for the right- 
whig test section. 


ground after the high-altitude portion of the flight is 
completed. The hot-film system, excluding the PCM, is 
estimated to weigh 20 lb and sits in the right whig under 
the hot-film sheet. 

Spectral analysis is the primary means of data 
reduction of the hot-film data. Preliminary calculations 
show that a 1-sec interval of data sampled at 20 to 
25 kHz is adequate to resolve the spectral content 
between 50 Hz and 1 0 kHz. This spectral content should 
be adequate for determining the flow field on the upper 
surface of the APEX airfoil. The detection of phase 
reversal and a significant change in power spectral 
density is expected to be the signature of laminar 
separation and the beginning of the separation bubble. 
Phase reversal of low-frequency spectra has been shown 
by Mangalam et al. to be an effective method of 
detecting laminar separation. 

Turbulent reattachment of the bubble is detected in 
the same manner as laminar separation — by phase 
reversal and a significant change in power spectral 
density caused by turbulence. The presence of vortex 
shedding is detected by performing both auto and cross 
power spectral density analysis on the hot films aft of 
the separation bubble. The detection of a significant 
increase in spectra in a specific frequency range and a 
consistent phase lag between the hot films is a signature 
of vortex shedding. 
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The Tolhnien-Schlichting instability wave frequencies 
in the free shear layer of the bubble is detected by 
spectral analysis of hot films inside the separation 
bubble and possibly by hot films just upstream of the 
separation bubble. The separated flow instability 
research of Dovgal et al. 4 show that instability waves, 
which cause transition, can be generated either upstream 
of the separation point or downstream of the separation 
point in the free shear layer. Their experiments show that 
harmonic disturbances causing transition exist inside the 
separation bubble. They discuss the concept of feedback 
interactions whereby instability waves are convected 
forward to the separation point as the flow circulates 
inside the bubble. 

At present, hot films have not been used to detect 
Tollmien-Schlichting instability frequencies and the 
bubble instabilities are assumed to be detected at the 
surface of the airfoil. The computational fluid dynamics 
analysis and Orr-Sommerfeld analysis performed by 
Tatineni and Zhong 8,9 suggest that the Tollmien- 
Schlichting instability waves occur at approximately 
1000 Hz. A significant increase in the spectral density in 
tins frequency range for hot films hi the separation 
bubble is, therefore, a measure of the Tollmien- 
Schlichting instability frequencies. 

Uncertainty Analysis 

Before actual flight, any estimate of the data quality 
from the hot-film system is difficult. The goal is a 
signal-to-noise ratio of 20 or greater. To reduce RFI and 
EMI noise, the hot-film strip comprises three laminated 
sheets. The top and bottom sheets are ground planes to 
shield the hot-film leads in the middle sheet. Twisted 
and shielded cabling are used for connections. The 
anemometry cards are packaged with ground plane 
protection. 

The aircraft power is filtered to ensure that the 
anemometry signals are not contaminated by power 
fluctuations. Special preflight ground test equipment is 
being developed that selectively blows both laminar and 
turbulent air over each hot film, matching the Nusselt 
number expected in flight. This equipment allows the 
individual hot-film signals to be compared and used to 
qualify, to first order, signal intensities between hot 
films. In addition, while the sailplane is suspended 
under the balloon, a 1-sec data interval is taken to assess 
noise levels. 


Vane. Total Temperature, and Accelerometer 
Measurements 

A pair of identical vanes is being developed for APEX 
to measure the angle of attack and angle of sideslip. The 
vanes will be calibrated to within ±0.25 deg to the local 
flow field angle of attack. The vanes will be mounted on 
a boom, located one-fuselage diameter in front of the 
nose of the aircraft. The position (upwash) error at this 
location, the difference between the local flow field 
angle of attack and the true angle of attack, can vary 
substantially. This difference in angle of attack is caused 
by the upwash from the boom, wings, and fuselage. The 
position error is expected to be a function of lift 
coefficient and Mach number. The error is estimated 
from file results of Rogallo 30 to be 0.3 deg at C t = 0.3 
and 1.2 deg at Cj = 1.2 for a Mach number of 0.6. A 
three-dimensional potential flow calculation of the 
APEX flow field provides a correction for the position 
error of angle of attack and angle of sideslip. 

The vane is currently being designed and fabricated 
and is expected to weigh approximately 10 gm to lower 
the moment of inertia and increase the vane response 
time. Preliminary near sea level wind-tunnel tests show 
that the vane damping ratios are between 0.2 to 0.33 and 
the natural frequencies are between 7 to 17.5 Hz as the 
dynamic pressure is varied between 18 and 50 KEAS. 
The in-flight time response of the vane is expected to be 
approximately 1.5 sec and is calculated by extrapolating 
the near sea level data to a flight condition of 40 KEAS 
at 100,000 ft by the method described in Barna and 
Crossman. 3 ’ The 1.5-sec tune response is adequate for 
steady-state measurements as flight simulations show 
that a flight condition of constant Mach number and 
constant lift coefficient for the sailplane can be 
maintained for 5 to 10 sec. 

The total temperature measurement is performed with 
an adaptation to a standard Rosemount (Burnsville, 
Minnesota) RTD total temperature probe. The RTD 
sensor is replaced with a Thermometries (Edison, 
New Jersey) thermistor to increase the time response of 
the sensor. Friehe and Kheltf developed this 
adaptation. The in-flight response time is expected to be 
approximately 1 sec and the accuracy is expected to be 
±0.5 °F. The sensor has been fabricated but has not been 
tested and calibrated. 

Three Endevco (San Juan Capistrano, California) 
piezoelectric accelerometers are mounted inside the 
wing to the upper surface at chord locations of 20, 50, 
and 80 percent. The accelerometers are sampled at 
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20,000 samples per second and have a frequency range 
of 1 to 10,000 Hz and an accuracy of ±0.05 g. A spectral 
analysis is performed on the accelerometer data and 
compared with the hot-film spectral analysis to explore 
any aeroelastic effects that may affect laminar-to- 
turbulent transition. 

Predicted Airfoil Performance 

The prediction of low-Reynolds-number airfoil 
performance is a formidable task that involves correctly 
modeling several flow phenomena as was shown in 
figure 2. Modeling the inviscid flow field including the 
presence of shock waves is generally considered the first 
step to determining the pressure distribution over the 
surface of the airfoil. The viscous flow field is composed 
of file boundary layer, laminar separation, laminar free 
shear layer, transition to turbulence in the free shear 
layer, turbulent free shear layer, reattachment of the 
turbulent free shear layer, and turbulent boundary layer, 
hi addition turbulent separation and laminar bubble 
separation, known as bubble bursting, are important 
physical characteristics to be modeled. The interaction 
between the inviscid and viscous flow fields can be 
significant. The presence of the separation bubble alters 
file effective shape of the inviscid airfoil. The classic 
assumption that pressure is constant across the boundary 
layer may not be valid across the separation bubble. In 
addition, boundary layers become large at low Reynolds 
numbers increasing the boundary displacement 
thickness, which can have an appreciable effect on the 
inviscid pressure distribution. 

The design and initial predictions of the APEX- 16 
airfoil are performed with the MSES code. The MSES 
airfoil design code uses the Euler equations to solve the 
inviscid flow field coupled with a two-equation 
dissipation integral method to solve for the viscous 
boundary layer. The transition location is determined via 

n 

the amplitude ratio ( e ) method, using growth rates that 
are precomputed from solutions of the Orr-Sommerfeld 
equation and correlated to the local shape factor 
parameter and momentum thickness Reynolds number. 
No compressibility corrections are employed partly 
because of the large uncertainty in the appropriate 
critical amplification parameter, w , for this high- 
subsonic-Mach-number and low-Reynolds-number 
flight regime. A value for n ■, of 12 was assumed for 
design of the APEX- 16 airfoil. Liebeck uses an earlier 
version of the MSES code, ISES, and finds that the 
predictions for low-Reynolds-number airfoils are 
reliable and accurate for low Mach numbers. One goal 
of the APEX experiment is to determine whether MSES 
remains reliable in the high-subsonic -Mach-number and 


low-Reynolds-number flight regime and what n jt 
values are appropriate. 

Figures 9 through 1 1 present the predicted drag polars 
and lift curves for the APEX- 16 airfoil for the chord 
Reynolds numbers of 200,000, 300,000, and 500,000, 
respectively, from the MSES code. The first apparent 
characteristic in the figures is the decrease in maximum 
lift coefficient with increasing Mach number. This 
decrease results from the separation of the turbulent 
boundary layer from the airfoil as the Mach number 
increases. The maximum lift coefficient decreases and 
the drag coefficient increases as the Reynolds numbers 
decrease. This result is expected as the separation 
bubbles become larger with lower Reynolds numbers, 
which decreases the overall performance of the airfoil. 
The lift curve slope is relatively unaffected by Mach 
number and Reynolds number except near stall. The 
slope of file pitching moment coefficients with angle of 
attack is also relatively unaffected by Mach number and 
Reynolds number. The predicted transition location, 
x TR , versus lift coefficient are also presented in the 
figures. The transition location on the upper surface 
moves forward and the transition location on the lower 
surface moves aft with increasing lift coefficient or 
angle of attack (a). 

Figure 12 presents the predicted drag polars and lift 
curves for various chord Reynolds numbers between 

200.000 and 2,000,000 for a 0.6 Mach number. The 
decreasing airfoil performance with lower Reynolds 
numbers is again apparent. An interesting feature of the 
figure is that the pitching moment increases with lower 
Reynolds numbers. Examining the data of McGhee 
et al., 34 the reverse would be expected. Their data show 
that, with lower Reynolds numbers, bubble reattachment 
occurs farther aft on the airfoil, which delays the 
pressure recovery on the upper surface. Tins effect 
results in decreasing the pitching moment. Figure 13 
shows the APEX- 16 airfoil predicted pressure 
distribution for a Reynolds number of 200,000 and 

300.000 at Mach 0.65. In the figure the bubble 
reattachment is predicted to move aft with lower 
Reynolds numbers and is in agreement with the data of 
McGhee et al. 34 The reattachment point is shown in file 
figure by the point of discontinuous change in slope of 
the pressure recovery on the upper surface. With lower 
Reynolds numbers, however, the separation point moves 
forward and the overall pressure on the upper surface 
increases. The point of separation is just after the point 
of minimum pressure on the upper surface as shown in 
the figure. The overall gam in upper surface pressure 
results in increasing the pitching moment with lower 
Reynolds numbers. 
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Figure 12. MSES prediction for APEX- 16 airfoil at M = 0.6. 
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Figure 13. MSES prediction of the pressure 
distribution over the APEX- 16 airfoil (M = 0.65; 
a = 3.5; Re = 200,000 and 300,000). 


predicted to become very erratic as the Mach number is 
increased into the transonic range, as shown in 
figures 14 and 15. The predicted interaction between 
tlie shock waves and the shedding vortices, as seen in 
tlie figures, has a profound effect on the flow field and 
tlie airfoil section lift coefficient. The section drag also 
increases substantially. A time-accurate Navier-Stokes 
analysis was also performed at NASA Dryden on the 
APEX-16 airfoil. Figure 16 shows the predicted 
unsteady separated vortex region on tlie aft upper 
surface of tlie airfoil. The results are similar to those of 
Tatineni and Zhong. 8, 9 


The vortex shedding criterion suggested by Pauley et 
al. 5 is 


P 


max 


®sep^ du\ 
u \dx) max 


-0.24 


( 9 ) 


A time-accurate Navier-Stokes analysis was 
performed on the APEX- 16 airfoil by Tatineni and 
Zhong. 8,9 Their analysis predicts that the separation 
bubble on tlie upper surface of the airfoil is unstable. 
The separation bubble is predicted to periodically shed 
at about 950 Hz for the subsonic flight condition of 
Mach 0.5, Reynolds number 200,000, and an angle of 
attack of 4 deg. The flow field over the upper surface is 


where P max is a dimensionless pressure gradient 
proposed by Gaster, 35 0 is tlie boundary layer 
momentum thickness at separation, n is tlie kinematic 
viscosity, and ( du/dx) max is tlie maximum velocity 


’Tatineni, Mahidhar and Xiaolin Zhong, “Numerical Simulation 
of Unsteady Low-Reynolds-Number Transonic Separated Flows Over 
the APEX Airfoil, APEX Critical Design Review,” NASA Dryden, 
1998, unpublished. Grant NCC 2-374, UCLA Flight Research Center. 








Figure 14. Unsteady variation of pressure contours for tlie transonic APEX-16 airfoil (M = 0.65; Re = 200,000; 
a = 4°). Time interval between frames is 0.0016 sec. 
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Figure 15. Unsteady variations of average section lift 
coefficient for the transonic APEX- 16 airfoil (M = 0.65; 
Re = 200,000; a = 4°). 


gradient. The shedding criterion for the APEX- 16 airfoil 
at a flight condition of Mach 0.65, Reynolds number 
200,000, and an angle of attack of 4 deg is 

P ma, = - 1-47 ( 10 ) 


Therefore, unstable shedding vortices should be 
expected on the APEX- 16 airfoil according to the 
suggested shedding criterion. 

The stability of the separation bubble has a large effect 
on the airfoil predicted performance. The MSES code, 
based on stable bubble calculations, predicts a lift 
coefficient of 0.96 at the flight condition of Mach 0.65, 
Reynolds number 200,000, and an angle of attack of 
4 deg. The Navier-Stokes code predicts an average 
section lift coefficient of 0.76 for the same flight 
condition. The Navier-Stokes analysis assumed laminar 
flow. The effects of turbulence on the stability of the 
separation bubble are uncertain. Gruber et al. 36 
performed a direct numerical simulation that showed 
that an amplified Tollmien-Schlichting wave in the free 
shear layer of a separation bubble develops into a large 
vortical structure. Whether the intensity of these vortical 
structures is large enough to maintain the structure as 
transition into turbulence occurs is unknown. These 
large vortical structures may be analogous to large-scale 
turbulent eddies that are quickly broken up in the 
turbulent flow field through vortex stretching and the 
tliree -dimensional effects of turbulent flow. 



Distance, in. Distance, in. 

’ ’ 990015 

Figure 16. Unsteady variation of velocity vectors for the aft section of the APEX upper surface. Plot sequence time 
interval is 0.0015 sec (M = 0.65; Re = 200,000; a = 4°). 
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Concluding Remarks 

The purpose of the APEX experiment is to increase 
the understanding of low-Reynolds-number airfoils in a 
low-turbulence flight environment. The APEX 
experiment regime is for altitudes between 70,000 and 
100,000 ft, Mach numbers between 0.5 and 0.65, and 
Reynolds numbers between 100,000 and 700,000. The 
following characteristics of the airfoil are to be 
determined: 

1. Section lift. 

2. Section drag. 

3. Location of the separation bubble. 

4. Vortex shedding characteristics. 

5. Tollmien-Schlichting frequencies. 
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